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THE PRICING OF LOOKBACK OPTIONS AND 
BINOMIAL APPROXIMATION 

KARL GROSSE-ERDMANN AND FABIEN HEUWELYCKX 
Universite de Mans 


Abstract. Refining a discrete model of Chenk and Vorst we obtain a 
closed formula for the price of a Enropean lookback option at any time 
between emission and matnrity. We derive an asymptotic expansion of 
the price as the number of periods tends to infinity, thereby solving a 
problem posed by Lin and Palmer. We prove, in particular, that the 
price in the discrete model tends to the price in the continuous Black- 
Scholes model. Our results are based on an asymptotic expansion of the 
binomial cumulative distribution function that improves several recent 
results in the literature. 


1. Introduction 

Cheuk and Vorst [ 3 ] have proposed a discrete model for the pricing of 
European lookback options with floating strike; they suppose implicitly that 
they evaluate the price at emission. In this paper we will refine their model 
in order to price the option at any given time after emission, and we derive 
an asymptotic expansion for the price as the number of time intervals tends 
to infinity. This solves completely a problem posed by Lin and Palmer [T3] ; 
in the special case of the price at emission the problem has already been 
treated by the second author [Ill- 

Lookback options give the holder the right to buy (for the call), respec¬ 
tively to sell (for the put), the underlying asset at maturity for its lowest, 
respectively highest, price during its lifetime. The payoff functions at ma¬ 
turity are therefore given by 

—min5i and max St —St, 
t<T t<T 

respectively. 

In the traditional continuous model (i.e., when the underlying asset price 
follows a Wiener process with drift, as proposed by Black and Scholes [2] and 
Merton [H]), Goldman, Sosin and Gatto m derived a formula for the price 
under the assumption that the spot rate r is non-zero. Babbs [T] obtained 
the price for r = 0 by passing to the limit. 

A discrete model for the price of a lookback option was proposed by Hull 
and White m, see also Hull |12j . who based it on the familiar binomial 
model of Cox, Ross and Rubinstein (CRR) [6], see Figure [LTl However, 
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since lookback options are path-dependent, Hull and White had to subdi¬ 
vide each node into different states. This problem was overcome by Cheuk 
and Vorst [3j who proposed an equivalent tree (CV) in which every node 
corresponds to a single state, see Figure [L^ (for a call). 



Figure 1.1. CRRtree 
with n = 3 



Figure 1.2. CV tree 
for the call with n = 3 


Neither Hull and White, nor Cheuk and Vorst obtained a closed formula 
for the prices in their model. Such a formula was first obtained by Follmer 
and Schied [S] and, by a different method, by the second author m (see also 
the discussion in m Appendix A]). However, all these papers only evaluate 
the price of the option at emission. 

The aim of this paper is two-fold. In Section 2 we will derive a closed 
formula, in the discrete model, for the price at any time after emission of a 
European lookback option with floating strike. In Section 4 we show that the 
price in the discrete model converges to the price in the continuous model, 
and we derive an asymptotic expansion. This answers a problem posed by 
Lin and Palmer m and generalizes the earlier results in m- In order 
to derive these asymptotics we need a refinement of the known asymptotic 
expansions of the binomial cumulative distribution function, which will be 
achieved in Section 3. The final section is devoted to numerical examples. 


2. Lookback options with floating strike 


In the discrete model of Cheuk and Vorst, the traditional CRR tree, see 
Figure fm is still used for modelling the evaluation of the underlying. 

We first introduce the usual notations for the various parameters: T is 
the time from emission to maturity, t with 0 < t < T is the present time, 
T = T — t is the remaining time until maturity, and St is the value of the 
underlying at time t. Moreover, r > 0 is the spot rate and a > 0 is the 
volatility of the underlying asset. 

Now, the binomial tree for the underlying is only built for the time interval 
[t, T] (the prices before t being known). This interval is divided equally into 
n subintervals. At each node, the price may either increase by a factor u or 
it may decrease by a factor d, where 

u = and d = u-^ = 


The probability of an increase is given by 


( 2 . 1 ) 


P = 


^rrln _ ^ 

u — d ’ 


we take n sufficiently large so that 0 < p < 1. 

Let us first consider t = 0, and let us concentrate on call options. For 
pricing a lookback option, Cheuk and Vorst have introduced a second tree 
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(CV), see Figure [L2l In the CV tree, the level j at time (0 < m < n) 
denotes the difference in powers of u between the value St^ at time tm and 
the lowest value of the underlying since emission. In other words, j is the 
non-negative integer such that 


As shown in [1] and made explicit in mi Theorem 2.1], the price of the 
option at emission is then given by 


where 

( 2 . 2 ) 


Cl‘{0) = So ^(1 - u->) Y, 

j=0 k=0 


q = = 


u - 

u — d 


(which lies in (0,1) if p does) and ^j^k,n is the number of paths in the CV 
tree from the initial node (0,0) to level j at maturity that have exactly 
k upward jumps. It was shown in m Theorem 2.1] that 


A 


j,k,n 


L-,) - L-U) 


else. 


In this paper we are interested in finding the price of a lookback (call) 
option at any given time t {0 < t < T). In that case, the market prices St*, 
0 < < t, of the underlying between emission and time t are known. Of 

course, the price at time t is not necessarily the minimal value; there is in 
fact an initial level jo that is defined by 

St = (min St *) u^° ■ 

^ t*<t ' 


In the sequel we will write for brevity 


Mt = min St* ■ 

t*<t 


Then we have 

(2.3) i„ = 

cry'T/n 

Clearly, Jq > 0 is not necessarily an integer, and the initial node is located at 
position (0,jo). This leads to a modified CV tree (see Figure [2Tl where m 
corresponds to time tm)- Note that after sufficiently many downward jumps 
one reaches a new minimum and thus level 0. From there on, all levels are 
integers. 

It can be shown as in [3] that 

n 

(2.4) Cl\t) = S, ^(1 - u-') E - 9)”"t 

jaJ k=o 

where J is the set of possible levels at maturity, q is given by (12.211 and, as 
before, A^.®^ ^ is the number of paths from the initial node (0,jo) to level j 
at maturity that have exactly k upward jumps. It remains to evaluate these 
numbers. 
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Figure 2.1. CV tree for the call with n = 5 and jo = 1-6 


6.6 


4.6 


3 

2.6 

2 

1 

0.6 

0 


Put abstractly, we have the following problem. Let jo > 0 be a positive 
real number and n € N. We create a graph with an initial node at (0,jo), 
and from each node {m,jm) at period m, 0 < m < n, we create two nodes 
at period m + 1 given by 

(m + l,jm + l) up 


and 


(m + 1, max(jm — 1,0)) down, 


with connecting edges. Let 

\jo 

j,k,n 

denote the number of paths from the initial node (0,jo) to the final node 
{n,j) that has exactly k upward jumps. The following result may then also 
be of independent interest. 


Lemma 2.1. Let jo > 0 and n G N. Then 


= 


n' 

kj ' 


n 

k- j 


n 


^ + LjoJ + 1 


n 

k-j-1 


if j = jo + 2A: - n 
and n — [joj < k < n, 

if j = jo + 2A; - n 
and n — < ^ < 

if0<j<n - [joJ - 1 
and j < k < [liiiiMiibij. 


All the other values are zero. 


Proof. (I) We will first consider the situation where jo is not an integer, see 
Figure 12.11 In that case there are three disjoint classes of terminal nodes. 
The proof will be done in three steps devoted to the nodes in each of these 
classes. 

(1) The first class is given by the end-points of simple solid lines. These 
nodes are exactly those for which the set of predecessors (going back to 
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m = 0) is the same as in a traditional binomial tree. It is then easy to 
count the number of paths. To arrive at such a node there are at least 
k = n — [joJ and at most k = n upward jumps, and there are (^) paths 
with exactly k upward jumps. The corresponding terminal levels are j = 
jo + k — {n — k) = jo + 2k — n. Thus, 


uo 



for j = jo + 2k — n and n — [joj < k < n. We note for later use that the 
levels j in this class range from jo + n — 2[joJ to jo + n in steps of 2, and 
that the number k of upward jumps is 


(2.5) 


k 


n + j - jo 
2 


(2) The second class is constituted of all the remaining terminal nodes 
that belong to the traditional binomial tree starting at (0,jo); in Figure [XT] 
they are indicated by dashed lines. We have here a partial binomial tree: 
some of the paths of the traditional binomial tree leading to such nodes have 
been lost. Indeed, some paths would have had at least one intermediate node 
at level zero. The possible terminal levels range from {jo} (if [joJ + n is 
even) or from {jq} + 1 (if [joJ + n is odd) to jo + n — 2[joJ — 2 in steps of 
two; moreover there are such nodes only when n > [joJ + 2. We can rewrite 
the lower bound in both cases as jo + n — 2[^^-^J. By (12.5p . each level 
corresponds to a unique number k of upward jumps, and the possible values 
of k for nodes of the second class range from n — J to n — [joJ — 1. 



Figure 2.2. The reflection principle 

Now, the number of paths leading to a node of the second class at level j 
is again but we have to withdraw all the lost paths. The number of 
lost paths is obtained by using the reflection principle, see Figure 12.21 It 
states that there is a one-to-one correspondence between the paths in a 
binomial tree connecting the nodes (0, a) and (n, b) (with a, b £ No) in 
touching or crossing the x-axis and the paths connecting the nodes (0, —a) 
and {n,b). Our case is equivalent to the situation where a = [joJ -|- 1 and 
b = j — {jo} -|- 1 = j — jo + LjoJ + 1- It remains to count the number of paths 
from (0, —a) to {n,b). The sum of the number of upward jumps {U) with 
the number of downward jumps (D) is the number of periods n, while the 
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difference between U and D is the overall increase, that is, a + h. We thus 
have that U = Therefore, the number of lost paths is 


n 

n+ffl+b 

2 


n 

n+2[jo\+j-jQ+2 

2 


n 


k + boj + 1 


where we have used (j2.5ll . Thus, 


\30 _ 


n 


+ LjoJ + 

for j = jo + 2k — n and n — < k < n — [joj — 1. 

(3) The third class of nodes is constituted by the remaining nodes; they 
are indicated in Figure ET] by double lines (either solid or dashdotted ones). 
They are the terminal nodes of a smaller Cheuk-Vorst tree with the initial 
node at (b’oj +1; 0). Moreover, there are such nodes only when n > [joJ +1. 
The possible terminal levels range from j = 0 to j = n — [joJ — 1. 

Unlike in the two previous cases, the paths joining the initial node (0, jo) 
with a terminal level j will not have a fixed number of upward jumps, just 
like in the traditional Cheuk-Vorst tree. We will prove that the number of 
paths with exactly k upward jumps arriving at a level j in the third class is 
given by 

' n \ ( n 

- jJ \k- j - 


\jo _ 


for 0 < j < n — [joj — 1 and j < A: < [ 


n-\jo\-l+j 


, and that no other values 


of k are possible. 

We prove the claim by induction on n, n > [joj -f 1. The result is trivial 
for n = [joj -|- 1. For n = [joj + 2, we can arrive at level j = 0 (with A: = 0) 
or at j = 1 (with A; = 1). In both cases, there is exactly one possible path. 

Now let n > [joj -|- 3. Again, for j = n — [joj — 1 (with k = n — [joj — 1) 
and j = n — [joj — 2 (with k = n — [joj — 2) the result is trivial. In both 
cases there is only one possible path, which are respectively [joj + 1 downs 
followed by n — [joj — 1 ups and [joj -f 2 downs followed by n — [joj — 2 ups. 

It remains to discuss the case when 0 < j < n — [joj — 3. If j = 0, there 
are two downward paths leading to it from period n — 1 if n -|- [joj is even, 
and three downward paths otherwise. 

(i) If n -|- [joj is even, they come from the nodes located at levels 0 or 1 
at period n — 1. 

If A; = 0, the only possible path comes from j = 0 and so 
\j0 _ \j0 _ 1 

^'■0,0,n ~ ^'■0,0,n-1 ~ 

If 1 < A: < then we have 

A io ^ A io I A io 

^^0,k,n ^^0,k,n—l ' ^^l,k,n—l 


n — 1 
k 
n 
k 


n — I 
k-1 


+ 


n — I 
k-1 


n — 1 
k-2 


n 

k-l 


n-[jo\-2 

2 


There are no paths for k > 
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(ii) If n + [joj is odd, there are additional downward paths coming from 
the partial binomial tree. By (12.5p . these paths did exactly 
ups before entering the Cheuk-Vorst tree. They make therefore no 
contribution if 0 < A: < and so we argue as in (i). 

If /c = it is impossible to come from level j = 0 at period 

n — 1 and so we have with (2) that 


A JO ^ A JO I A JO 

^0,/c,n ^l,k,n—l ^ {jQ},k,n—l 

'n-1 
^k-i 
''n-l 

''n\ / n 


A:- 1 


n — 1\ /n — 1 

A:-2y ^ V k 

n — 1\ fn — 1 
k-2)^\ k 


n — 1 

k + [jo\ + 1 

n — 1 
k-1 


There are no paths for k > ^ 

If 1 < j < 77 ~ LjoJ ~ there are two paths leading to it from period 
n — 1, one upward path coming from the node located at level j — 1 and one 
downward path coming from the node located at level j + 1. If 0 < A: < j 
there are no paths with k ups. If A: = the only possible path comes from 
j — 1 and so 


A JO ^ A JO 


= 1 . 


If J + 1 < A: < |^!L_LZ2J_Ji:ij ^ then we have 


A JO ^ A JO , A JO 


n — 1 
k-j 
n 

k-j 


n — 1 \ ( n—1 

+ 


k — j — ij \k — j — 1 
n 

k-j -I 



There are no paths for k > , 

This proves the lemma in the case when j^ is not an integer. 

(II) The situation changes slightly when jo is an integer, see Figure [T3l 
The three classes of nodes are defined as before: Those for which the set 
of predecessors is the same as in the traditional binomial tree starting from 
(0)io) (indicated by solid lines); the remaining terminal nodes of the tradi¬ 
tional binomial tree (indicated by dashed lines or by thick solid lines); and 
the terminal nodes of the smaller Cheuk-Vorst tree starting at (jo -|- 1,0) 
(indicated by double lines, thick solid lines, or the dash-dotted line). Note, 
however, that some terminal nodes belong to both the second and the third 
class. 

(1) For the first class, we apply the same argument as before. 

(2) For the nodes in the second class we only count the paths that are 

inside the binomial tree. As before, there are such nodes only when n > 
jo -|- 2, the possible number k of upward paths ranges from n — to 

re — jo — 1, and the terminal level is given by j = jo -|- 2A; — re. The number of 
corresponding paths (inside the binomial tree) is (^) minus the number of 
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Figure 2.3. CV tree for the call with n = 5 and jo = 1 


lost paths. But now the lost paths are those paths in the binomial tree that 
hit level —1 at some point. This is equivalent to counting the paths from 
(0)Jo + 1) to (n,j + 1) that hit level 0, and the number of such paths is, as 
before, • This confirms the lemma for the stated values of j and k. 

(3) Finally, for the nodes of the third class we count those paths that at 
least once follow a path outside the binomial tree (in this way we do not 
count paths twice when the node also belongs to the second class). The 
argument as before confirms the third alternative in the lemma. 

We note that if the terminal level j belongs to both the second and the 
third class then the paths inside the binomial tree have exactly k = 
upward jumps (whence n + j — jo is even), while the paths outside this tree 
have at most 


n- jo -1 + j 


< 


n + j - Jo 


upward jumps. This shows that for such levels j there is no conflict in the 
statement of the lemma. □ 


Combining the formula (12.411 for the price of a lookback call with Lemma 
O we obtain the following. 

Theorem 2.2. Let 0 <t <T and n € N. The price of a European lookback 
call option with floating strike at time t is given by 

( 2 . 6 ) cl\t) = St{V,-V2 + V3), 

where 


U—U \^/ 


k=k^ 
^-LjoJ-1 


h—h ■ 

— '^min 

fl-LioJ-l fcmax 


V2= ^ (1- , )g^(l - gr~^ 


k + [jo\ + 1 


j=0 


k=j 


n 

k-j 


n 

k- j -I 


q\l-q) 


n—k 
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with jo given by ()2.3p . q given by (12.2p . fcmin = n — and femax = 



Figure 2.4. CV tree for the put with n = 3 

In the same way one can determine the price of a lookback put option 
with floating strike. In that case we would have negative values for levels 
because the underlying price is always lower than the maximal price, see 
Figure [231 For convenience we will write these levels as —j with j > 0. In 
terms of the levels j we are then led to the same tree as in the case of the 
call. 

With this adaptation, the level for the put at time tm is given by a non¬ 
negative number j such that 

St^ = ( max St*)u~^. 

The initial level jo (at initial time t) of the tree satisfies 

St = ( max St *) u ~^°, 

^ t*<t ' 


so that 

log(maxt*<f St* I St) 
a*/^r/n 

It can be shown as in Cheuk and Vorst [1] that the price at time t of the 
European lookback put option with floating strike is given by 

n 

= StY.in’ -i)Y. - vt t"-\ 

jaJ k=o 

where J is the set of possible levels (taken positively) at maturity, q is 
given by (j2.2p and ^ is the number of paths from the initial node (0, jo) 
to level j at maturity that have exactly k upward jumps. Since we have 
evaluated these numbers of paths in Lemma l2.II we obtain the following. 


Theorem 2.3. Let 0 < t <T and n G N. The price of a European lookback 
put option with floating strike at time t is given by 

Pl\t) = SflVi -V2 + F 3 ), 
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where 


= X] 

h—h ■ 

— '^min 


jo+2k—n 


n 


-i)( 


k ^n—k 


.jo+2k—n 


- 1 ) 


V2= ^ (U- 

h—k • 

— '^min 

^ LJoJ 1 ^max 

V4 = (rf - 1) E [ 


n 


^ + LjoJ +1 


{l-qfqn-k^ 


n 


n 


k-jj \ k - j-l 




k ^n—k 


j=0 k=j 

with jo given by (|2.7p . q given by (12.2p . fcmin = n — and /cmax = 


3. Asymptotic expansions of the binomial cumulative 

DISTRIBUTION FUNCTION 

In order to obtain asymptotic expansions for the price of lookback op¬ 
tions as the number of steps n tends to infinity we first need to derive an 
asymptotic expansion of the binomial cumulative distribution function 

jn / 

^n,pn{jn) = 

k=0 ^ 

or, equivalently, for the complementary binomial cumulative distribution 
function 

n , 

^n,pSjn) = 

k=jn + l ^ 

and that with a lower error term than in the expansions found in the litera¬ 
ture. Moreover, our method allows us to treat any sequence (pn) satisfying 
the rather weak assumption that 

0 < liminfpn < limsupp„ < 1. 

n—>-00 

Our approach relies on the work of Uspensky [16] who represented bino¬ 
mial probabilities in a convenient analytical form. Chang and Palmer |3l 
Lemma 1] refined this result in order to obtain an approximation of order 
o{n~^) in the case when pn 1/2. Recently, Lin and Palmer [141 Lemma 
C.l] provided an estimate of order 0{n~‘^) if = 1/2 -|- 

Chang and Palmer have used their approximation in order to obtain 
an asymptotic expansion for digital and European options with a remainder 
term (see m Theorem 1.1]), while Lin and Palmer [T^j treated 

barrier options with the same precision. The second author m recently 
evaluated lookback options with floating strike, limiting himself to the price 
at emission (that is, at t = 0). Moreover, when the spot rate r equals 0 then 
he only obtained a remainder term 0{n~^) (see [111 Remark 4.1]). 

In order to obtain an asymptotic expansion for lookback options with 
floating strike with a remainder term 0(n“^/^), and that for any spot rate 
r > 0 and for the price at any time t > 0, we will need an approximation 


PnC^-Pn) 


n—k 


)pn(l - Pn) 
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of the binomial cumulative distribution with remainder term as 

provided by the following theorem. 

In the sequel, C > 0 denotes a generic constant, which may have a differ¬ 
ent value at each occurrence. Note also that, for the sake of readability, we 
will often drop the index n. 

Let $ denote the standard normal cumulative distribution function. 


Theorem 3.1. Suppose that p = pn satisfies 


0 < liminfp^ < limsuppn < 1. 


n^oo 


n^oo 


If q = qn = '^-Pn and 0 < j = jn< n, then 


E 


fc =0 


'a\ u ri — k 

pkqn k 


= $(2/) + 




e 2 


ft 


/ ft ft 
\y/v F3/2 


+ ^) +0 


n5/2 


as n ^ oo, where V = npq, y = ^ and 

ft = l{q-p){^-y‘^), 

ft = 2/[^(-3 + 7y2 - /) - i(-3 + lly2 _ 2y% 

P3 = {q- ft [^(123 + 129y2 _ 384/ + 95/ - 5/) 

- 3io(3 + 69/ - 399/ + 145/ - lO/)], 
ft = y[Tdbo(-4293 - 1359y2 + 6165/ - 1971/ + 185/ - 5y^°) 
+ 3^(3106 1395y2 _ 7794/ -h 2979/ - 325/ 

+ ^fe(135 - 1035/ -h 7947/ - 4167/ -h 560/ - 20y^O)]. 


As mentioned above, the proof of this result will be based on an analytical 
representation of binomial probabilities due to Uspensky m p- 121]. 


Theorem (Uspensky). Let 0<p<l, q = l—p and 0 < j < n be fixed 
numbers. Then 


where 


y = 


j -np+\ 

/U 


and y' = 


np+\ 


with V = npq; here the function J is defined by 


(3.1) 

where 


Ay) 


J_ r n sm{yVV(f-x) 
2-71 Jo ^ sin f 


dyj. 


y G R, 


p = + q\, oj = arg(pe*‘^ + q) and x = 'aoj — npp. 

We will derive Theorem 13. II from Uspensky’s representation by employing 
and refining his ideas (see m pp. 121-129]). For this we will need two 
preliminary lemmas. 
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Lemma 3.2. Suppose that p = Pn satisfies 


0 < liminfp^ < limsuppn < 1- 

n—>-oo 72 —)-oo 

Let q = Qn = I — Pn and V = npq. For a fixed constant M > 0, let ip be a 
positive number such that p < . Then we have for R = + q\^ 

that 

R = (l + R^p* + Rep^ + ^Rjp^ + 0(nV^^+®)) 

k=l 

as n ^ oo, where 

R2 = -^V, 

^4 = \V[l -pq), 

- iP^ + pV)- 

Proof. Let 

p= Ipe'-'^ + q\, 

so that R = p^. We then have that 

p = {p^ + ‘^pq cos p + q^Y^‘^ = — Apq sin^ ^ . 

This gives us that 


(3.2) 


1 / 2 
log /^ = X log 1 “ sin - 


= -2pqsiT? I - A{pqf sin^ | - ^{pqf sin^^ f “ 


where 5 = 32(1 — r/)“^(pg)^ sin® ^ for some real number p between 0 and 
4pgsin^ Since liminf„^ooPn) hminf^^oo > 0 we have that V = npq —)■ 
oo and therefore 0 < (^ < —)■ 0 as n ^ oo. Thus 

(3.3) 0 < (5 < Cp'^q'^p^ 

for n sufficiently large. 

Now, in order to obtain bounds for log p, we use suitable Taylor expansions 
about 0. First, we have 


■ 2 T 1-2 1- 4 , 1- 

sm — = —p - p -\ - 

2 4^ 48^ 1440 


Thus we have 

(3.4) 4^ J_ 

^ ’ 4^ 48^ 1440 


80640 




1 


7257600 


p^^ + o{p^^). 




- p < sin — < —p - p 4- p 

80640^ “ 2 “ 4^ 48^ 1440^ 


for n sufficiently large. Next we deduce from 


■ 47^ 1-4 1-6| 1- 8 

2 16^ 96^ 1280^ 


17 


p^^ + o{p^'^) 


lOl 


that 

(3.5) 


483840 
1 


1 


— p^ - p^ < sin^ — < — p"^ - p'^ 4- 

16 ^ 96 ^ “ 2 “ 16 ^ 96 ^ 1280 




for n sufficiently large. Finally, we obtain from 


1 


1 


7 


sin — = —p - p 4- 

2 64^ 256^ 15360 


p^^ + oip^'^) 


lOl 
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that 

(3.6) 


1^6 _ ^ 8 < V < 1 6 

64^ 256^ “ 2 - 64^ 


for n sufficiently large. From 
lower bounds for the expression 


3D, (|3.5p , (13.6p . we hnd upper and 


1 


1 


logp<-2pg[-^ -48^+1440 




80640 




(3.7) 


and 


(3.8) 


7 I 4 

7 6\ 

3 3 / 

' 7 6 

7 8\ 

(le^^ ■ 


- yP q ( 

■ 



-^2 2 , 7?4 4 , TZe 6 , 7 

= —p + —if + —^ + —pq 
n n n 24 


1680 




log p > —2pq( — p^ -H-y?®') — 4p^q^ (—- p^ 

sp- 48^ 1440^ / ^^Vl6^ 96^ 


+ 


1 


1280' 


32 




1 


.,4 4, 8 


R2 


I 2 , R'i A , Rq 6 22/7 |,^2 2\8 

-p + —p + —p -p q [—+Cpq ]p . 
n n n V320 / 

By combining (13.7h and (|3.8D . we establish that 

gR2¥5^+R4VP'‘+R6¥5®gAl ^ _ pTl ^ gR2‘P^+R4¥’'*+J?6‘P®gA2 

where 

= -Mwo +^ = YAiL^’’'^y 4 0. 

This clearly implies that 

1^ _ gR2‘^^+^4V5"‘+H6<P® I < gR2‘P^+R4‘P'*+R6¥’® ^gA2 _ 

Using the facts that < 1 + 2x for 0 < x < 1 and e* > 1 + x, we have for 
large n that 

gA2 _ gAi < 2 A 2 - Ai < CVp^- 
note that A 2 —)• 0 as n ^ 00 . It follows that 

1^ _ gR2¥’^+R4</’'‘+R6¥’® I <; (jQR2<P^+R4‘f'’‘+R6V^Yp^ 
in other words, 

(3 9) jl — gR2v‘^+R4v^+R6v^ _|_ Q^^R2‘P^+R4v‘^+R6‘p^Y^Sy 

Now, the definitions of i ?4 and Rq and the fact that 0 < imply 

that 

(3.10) R4P^ + Rgp^<C. 

By combining (|3.9p and (|3.10l) . we see that 

(3.11) R = +0{np^)). 

We rewrite the second exponential term with a Taylor expansion about 0 
to obtain for some p between 0 and R^p^ + RqP^, 

^R4‘P +R6¥’ = 1 -|- (i?4</?^ + R&P^) + 7t(7?4<7’^ + + —A{R4p'^ + RqP^)^, 

2 6 
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hence with (j3.10l) 

(3.12) = 1 + i?4/ + + 0(nV^°) + 0(nV^^)- 

Combining (I3.1ip and (|3.12p we obtain the result. □ 


Lemma 3.3. Suppose that p = Pn satisfies 


0 < liminfpn < limsupp^i < 1. 


n—)-CxD 


n^OD 


Let q = Qn = I — Pn and V = npq. Then, for any non-negative integer m, 
there exists a constant C > 0 such that 

/•OO 

1^:= / 

Jo 


Moreover, suppose that r = > 0 satisfies r ^ = 0{n°‘) with a < 1/2. 

Then for any integers m and k, there exists a constant C > 0 such that 


/ OO 


Proof. We first note that by our hypotheses there is some (^ > 0 such that, 
for large n, 

(3.13) V > dn and tW > 

We set X = y/Vip. Then the integrals become 

poo 

/^ = l/-(”^+b/2 / e-2^\'^dx 

Jo 

and 

poo 

/;; = y-(”^+b/2 / e-2^\^dx. 

J Ty/V 

The integrands can be bounded by a function Ce~^ if m > 0, or if m < 0 
and X > 1. Since, by (I3.13p . ry/V —OO as n —OO, we obtain that 

lyn < 

and 

I*m < 

Thus we deduce with (j3.13p that 

Im < 


and, for any k, 

Rm < Cn-\ 

□ 


We can now prove our main result. 

Proof of Theorem \3.1\ The proof will be divided into several steps. 

(1) In view of Uspensky’s representation stated above the value J[y) plays 
a crucial role, see m- We will split the integral into two parts. Let 

r = 
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SO that r ^ Throughout the proof we suppose that n is large 

enough to have that t < f • Then 

1 


r{y) ■■= 


2tt 


i: 


sinf 


< 


I r p- 

‘T / 

l-K Sin 


With (j3.2p we have that 


9 (/? 

n log p < —2V sin —. 


Applying the fact that sin % > — for 0 < (/? < vr, we obtain that 


J 


1 f 


TT g-2y sin2 I 


sin ^ 


dip < 


1 r 
^Jr 


2V,.2 1 




< K 
ip 2 


2V,.2 1 






dp. 


Substituting x = -p and applying Lemma [331 we have that 


(3.14) 


"'W S 2 L 


i 1 

e-iV^-dx = 0 

X 


5/2 


n 


since (pr) ^ From ()3.ip and ()3.14l) we then obtain that 

(3.15) J(,) = ±/'T..=(A%^d^ + 0( ' ). 

271 Jq sm ^ 

(2) Looking at the integrand of J{y) we now want to estimate sin(o — x) 
in powers of p, where a G M; recall that 

uj = argPpe*"^ + q) and y = nu — npp. 

By Taylor expansion we have that 


(3.16) 


sin(a -x) = sm(a) - cos(a)x - ^ sin(a)x^ + ^ cos(a)x^ 

2 0 

+ ^ sin(a)x'^ - — cos(a - p)x^ 


for some p between 0 and x- Since by (I3.15P it suffices to assume that 
0<p<T<^ we see that 


psin<y9 

UJ = arctan-, 

p cos p + q 

so that 

psinp 

X = n arctan - npp. 

p cos p + q 

Now, 

d / x(/g) \ 1 p-q _ I _ 

dp\ n ) 2 ^ 21 + 2pq{cos + — 1) ’ 

where we have used that p + q = 1. Note that, as x ^ 0, 


——- = 1 - 2pqx + {2pqfx‘^ + 0{x^), 

1 + 2pqx 

where the constant in the big-0 condition does not depend on n since pq is 
bounded in n. Thus the Taylor expansion 
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gives US that 

d 

dip\ n J 


1 , P-<1 

= 2 




l-2p9(-Y + fj)+(2P'!t^+0(4.' 


= ^pq{p - - :^pq{p - 9)(i - + o{ip^) 

X = X3<p^ + + nO(v9^) 


^pq{p - q)'P^ - ^ 

and hence 
(3.17) 
with 

X3 = \y{p- q), 

X5 = -m^(p ~ 

Applying (I3.16p and (|3.17l) we obtain 
sin(a-x) = sin(a)-cos(a)(x3<^^+X57^^+nO(v7^)) 

- X sin(a) (xsV^+X 5 (p^+nO{(p'^)f + ^ cos(a) {x 3 (p^+nO{(p^)f 
2 0 

+ ^sin(a)(x3V9^+nO(<y9^))‘^-^cos(o-r/)(nO(v9^))^ 


and hence 


(3.18) 


sin(a-x) = sin(a)-cos(a)x3<7’^-cos(a)x57^®-2 sin(a)x3V^^ 
-sin(a)x3X57’^ + ^ cos(a)xi¥?® + ^ sin(a)x37’^^ 


+ ^0(nV^^+^)- 


k=l 

(3) Next, using Laurent expansion, we have that 

12 1 7 

-£ ~ T + “I- 


sin f ip 12' 


2880 


+ 0{p^), 


which together with Lemma 13.21 and the fact that R = gives that 


= 


Sin • 


(3.19) 


2 1 
p 12 


7 


+ 2ii4 


-2880 


7^^ + 


1 

L12 


i?4 + 2i?g 


<P 


k=0 


We can now rewrite the integrand of (|3.15p . Setting 

yVV p = ap = a 

we obtain, by combining (I3.18P and (j3.19l) . 

„ sm{yVV p-x) 


(3.20) 


sinf 




11 


sm a 


7^) + ^ 'hp^ + ^ 


k=l 


k=l 


where Jio = 0 and 
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Ji = ^ sin(cK^), 

J 2 = -2x3Cos(a(^), 

•^3 = + 2 ^ 4 ] sin(av9), 

’h = - [i^X3 + 2 x 5 ] cos(a(/7), 

J5 = [^-^4 + 2/^6 - xi] sin(Q;(y9), 


Je = -2i?4X3Cos(a(^), 

J 7 = [Rl - 2x3X5 - ^xi] sin(a 99 ), 
-^8 = 5X3COs(q;(/?), 

Jg = -i?4xisin(Q;(y9), 

•^11 = n>^3sin(a(^). 


(4) We will estimate (13.151) using the form of the integrand what we have 
obtained in (I3.20p . 

We begin with the error terms. Applying Lemma 13.31 we have that 

(3.21) I £ d<^| < Cn’^ dy. = O (^) , 

with fk = k an integer between 1 and 8. 

For the main part in (I3.20p we can replace the integral on [0, r] by one on 
[0,oo) with an error of at most 0(n“^/^). Indeed, 

/ OO / o . 

gin(Q;(/?) + ^ 

^ k=i 

11 




Using the definitions of the coefficients and noting that = 0(n^/^), 
we obtain with Lemma 13.31 that 

(3.22) r*{y) = o{-^ 

Vi 


or 


5 / 2 ;- 


Applying ()3.20l) . ()3.21l) and (|3.22l) to (I3.15p . we derive that 

1 /'OO o 1 31 /.Qo / 1 \ 

■'<!'> = 27/0 ^ ^ + 0(7572 ) 

k=l 


with a = ?/\/U and R 2 = 

It remains to evaluate these integrals, which we have relegated to the 
Appendix® After simplification, using in particular that {p — qY = 
we obtain that 


(3.23) J{y) 


1 e 2 ?^ 




/ 'Pi(y) 'P 2 (y) P^{y) PA{y) \ I p/ M 

V Vu 1^ ^ U3/2 u2 Vn5/2; 


with 

'Pi(y) = g(g-p)(l -y 2 ), 

^ 2 ( 1 /) = 2/[^(-3 + 7y2 - - 2|(-3 + lly2 _ 2y% 

PM = (g-p)[6^(123 + 129y2 - 384/ + 95/-5/) 

- 3io(3 + 6V - 399y4 + I45y6 _ l0y% 

My) = y[ 15 ^(-4293 - 13597/2 + 6165/ - 1971/ + 185/ - 5/°) 
+ 3^(3105 + 1395/ - 7794/ + 2979y® - 325/ + 10y^°) 
+ ^^(135 - 1035/ + 7947/ - 4167/ + 560/ - 207/1°)]. 
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(5) Finally, by Uspensky, we have that 


(3.24) 

where 


k=0 




j -np+ I , np+\ 

y = --=—^ and y = -- 


Vv ^ Vv ' 

Since 0 < lim inf^^oo Pn < 1™ < 1 there are J > 0 and C > 0 

snch that, for large n, 

6n < V < n and 6^/n < \y'\ < C^/n. 

It follows that for each k = 1,..., 4 


(3.25) 


e Pfc(y') /I 


n 


5/2 ) ■ 


Moreover, for x < —2, — < x. Therefore, if n is sufficiently large, then 


(3.26) 

so that 

(3.27) 




Now the theorem follows from (j3.24|] . (13.2311 and (I3.27p . 
We can deduce a first easy corollary. 


□ 


Corollary 3.4. Suppose that p = Pn ^ Po with 0 < po < 1 o,nd that j = jn 
satisfies ^ jo- Then 


t (:)/(!-.)”--0(4 


k=0 

j 


5/2 


k=0 




n5l2 


if jo < Po, 
if jo > Po- 


Proof. From our assumptions it follows that 


y 


j -np+ ^ 


30-Po 


Pn ny/p{l -p) v^po(l -po) 
Then as in (|3.25|) we have that, for A: = 1,..., 4, 

1 




= O 


n 


5/2 


^ U ^/2 

Moreover, if jo < po then as in (13.261) we have that 

1 

on® 

while if jo > Po, we have that 




1 
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The claim now follows from Theorem ixn □ 

The second corollary will be crucial in our applications to lookback options 
in the next section. 


Corollary 3.5. Suppose that 


1 


a 


+ '1 + 


7 


H—^ + 


^ 2 ^/n ' n ' n3/2 ' 


n 


n5l2 


+ 0 


ri'- 


and 


where {bn)n is a bounded sequence. Then 


j-- + aVn + - + 6„ + ^ + - + ^ 


+ 0 ^ 




E 

k=j 




Bn ^ C 0 -C 2 BI 


n 


Do — DiBn — D^Bn Eq — EiBn — -E' 2 -Bn + D 4 B. 

' o /o ~r o 


n 


3/2 


4 

n\ 




n 


5/2/’ 


where 


Co = 2a^A - (1 - A^)(A - 8a)/12 + C, 

C 2 = A/2, 

Do = AaPA + 2(1 - A^)p/3 + D, 

Di = {8aA - l )/6 - (1 - ^ 2)(^2 _ 24a2 _ 3)/i2 + AC, 

Ds = {l- A‘^)/6, 

Eo = 2(/32 + 207 )^ + (1 - A^)(iba^C + 2^)/3 + (3 - A^){iba^ - 2C)aA/3 + 
- 4yl2 + 1) (iGa^ - C')/12 - (5A® - bSA"^ + 33^2 + 171)21/1440 + 
(52l6-4l2l^+2l2l2+27)a/90-(724^-402l2 + 15)a2^/18-2lC'2/2+.E, 

El = AI3A/3 + (1 - ^2)(27l _ i2a)/3/3 + AD, 

E 2 = 2a2^+(i_yl2)((7_^2a2yl)/2-(^4_g^2_^9)^/24+(^^-6A2_^3^Q,/3^ 

E4 = {3- 712)71/24, 

with A = 2{a — a), Bn = 2{j3 — bn), C = 2(7 — 0 ), D = 2{5 — d), E = 2(e —e). 


Proof. The proof is based on the same ideas as Lemma 3.1 of m- As it 
is very computational, we will just recall the main steps for obtaining the 
result. 

We first note that by Theorem 13.11 


E 

k=j 


i -1 

p"(l-p)"-" = l-E 

k=0 


n 


p’^ii-p) 


n—k 


= 1 - $(y) - 






Pi P 2 
—- + — + 


^ 3/2 ^ Y2j ^^y.^5/2j 


with y = ^ ) where Pi,..., P 4 and V are as in that theorem. The 

goal is then to give an asymptotic expansion of each term in this sum. 

We begin with 1 — <h(y) = <h(—y). This term is an integral that we will 
decompose into two parts: one from —00 to A (giving ‘h(A)) and second one 
from A to —y. This splitting is motivated by the convergence of y to —A 
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as n tends to infinity. Then we find an asymptotic expansion for the second 
integral by using a Taylor expansion about A in which we substitute y by 
its asymptotic expansion. 

We finish with all the other terms. For each of them, we substitute V 
and y by their respective asymptotic expansions. □ 


Remark 3.6. We have stated the result for 


^n,pU - 1 ) = E 

k=j 


pHi-pY 


in line with the results in and m- The corresponding result for the 

binomial cumulative distribution function 

k=0 ^ ^ 


is obtained by writing j in the form 

n 1 c d e 

and taking 1 minus the result given in the corollary. 



4. Asymptotics of the price for lookback options 


In this section we combine the results of the previous sections in order to 
derive asymptotic expansions for the price of lookback options. We will use 
the notation of Section 2. In addition we adopt the following notations that 
are in line with common usage in the literature: 

di = ^(log^ + (r+4)^)> 

d2 = ^(log^ + (r- 4)^) =di-a,/^, 

ds = ^(-log^ + (r-4)^) =-di + ‘^V^. 

^4 = ^( - log ^ + (r + ^)t) = ds + 
and we will write 
Kn = {jo}(l - {jo}), 
where, as before, jo = 

r/n 

We first note that the price Cn of a lookback call as a function of n shows 
some (mild) oscillations, see Figure 14.11 which are caused by the fact that 
the non-integer part {jo} of the initial level jo varies with n. 

This reminds one of the (much more violent) oscillations for European 
vanilla options. To deal with them, Diener and Diener [7] have used asymp¬ 
totic expansions of the form 


a = co + ^ + - + 0 
'n n 


n3/2 


where the coefficients ci et C 2 are allowed to be bounded functions of n. The 
variability of these coefficients capture the observed oscillations. 

In order to obtain such asymptotic expansions in our setting we will need 
the following. 
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26.4 

26.3 

26.2 

26.1 

26 

0 100 200 300 400 ’ 60 80 100 120 140 

Figure 4.1. Price of a lookback call, and its fine structure 



Lemma 4.1. Let oq, ui, 02 £ o-nd let rj = ri(n) be a bounded function 
of n. Then 


1 + 


ai 


+ 


n 


a2 


n 




= 1 + 


aiv «2?? - |??(1 - i?)Qi ^ 
'n n 




Proof. We have that, for large n. 


l + ^ + “2 

n n 




. 71^/2 


= exp 

= exp 

= 1 + 


01 ^ 02 ^ ^ 


aiTj 


n n 

(02 - ^al)r] 


1 


+ 


+ 0 


n 


n 


n ^/2 
1 


n 


3/2 


aih ^ (02 - \al)ri ^ Ic^ ^ ^ 
fn n 2 n 


n3/2 


which confirms the assertion. It is important to note here that each big-0 
condition contains a constant that is absolute. □ 


In essence, the lemma says that we may expand the function (1 + + 

^ as if rj was a constant. Diener and Diener therefore speak 

of a frozen parameter. 

We first consider the asymptotic expansion of the lookback call. Its price 
in the Black-Scholes model is well known. If r > 0 then Goldman, Sosin and 
Gatto m found that 

(4.1) Ci^^it) = St- StOiB^ - MtB2 + St{l- 02)^3, 

where 
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B2 = e-""$(d2), 

^3 = e— 

By passing to the limit in (14.11) . Babbs [T] obtained the price in the case 
r = 0 as 

= St- StBi - MtB2 - St{B*^ - Bl), 

where 

u* _ ^ n: 

-84 - O-yV ^ . 

Theorem 4.2. Let 0 < t < T and n G N. The price at time t of the Euro¬ 
pean lookback call option with floating strike in the n-period CRR binomial 
model satisfies the following: 

(i) if r > 0 then 

1 


= CfJsit) - St^{e,B, + 92 Bs) 


Si 


+ 0 


a‘^T 

IT 

1 


2 ' ' y/n 

01 + 2)Bi + (02 + 2 - TOBa) 


— MtT2Bi 


1 

J n 


n 


•il2 


where B 4 = ^ ^ ^T T = ^ 02 /tn - (1 + ^) log ^ 

and r 2 = i + log %-■ 


2 ' ' 6cr77 

(ii) if r = 0, then 


CL\t) = - St^{2B, + Bl - Bl)^ 


a‘^T 




St— ((3 + 3tt„ - —)Bi + Bl) - StTfB^ 


J n 




where T| = ^ + + 


12 


i^log A 


6 (TyG- 


Mt- 


Proof. The hrst step of the proof consists in writing (12.6p as a combination 
of (complementary) binomial cumulative distribution functions. 

As for Vi we note that, by (|2.1I1 and (12.211 . 

^n-jo-2k^k _ ^^n-k ^ _ ^^n-k^ 


SO that 


Fi = b;,(ji - 1) - ^e—B;,(ii - 1), 


where ji = n — . 

For F 2 we first make a change of index k ^ k-\- [joJ + 1 and then proceed 
as for Fi to obtain that 

F 2 = q-l^°J-'b;,(j 2 -1) - ^e-’'"p-L^°J-'B;p(j2 -1), 
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where 


and j 2 = ji + boj + 1 - 

For V 3 we proceed as in m- We first split the inner sum and then 
change indices, k ^ k — j and k k — j — 1. Next we interchange the 
two double sums that have appeared, noting that 0 < j < n — [joj — 1 and 
0 < k < [Hzizlioinij is equivalent to 0 < A: < J and 0 < j < 

n — boj — 1 — 2A:; and that 0 < j < n — [joJ “ 1 and 0 < k < [IlzizlMnij _ 1 
is equivalent to 0 < A: < J _ i and 0 < j < n — boJ — 3 — 2A:. Note 

that = ji-1. 

Now, if r > 0, then the geometric series that arise as inner sums have 
ratio different from 1 , so that by continuing as in m we obtain 

^3 =(Q ^i)(Qd-l) - QBnAh - 1 )) 

n-boj-i 

+ (Q^n,l-,(J3) - i3„,l-,b3 - 1 )) 

lOdWLjol-l 

+ b - Qd) 

with j 3 = ji — 1, where we have used that u~^ = d. 

However, if r = 0, then Q = uso that geometric series with ratio 1 appear. 
Using A:(^) = the calculations then lead us to 

^3 =(boJ -n - {^n,q{h) - uBn,q{h “ 1)) “ ‘^uBn,q{h “ 1) 

^-lioJ -1 

+ ^ {uBn,p{h) - ^n,p{h - 1)) 

+ 2nq{Bn-l,q{h - 1) - uBn-l,q{j‘i “ 2)), 

where we have used that 1 — g = p in this situation. We note that the 
formulas for Vi and V 2 are not affected, but one may replace Q hy u and P 
by d in V 2 . 

The second step of the proof is to expand each term using Corollary 13.51 
and Remark 13.61 We note that we use Lemma O with 7] = {jo} to expand 
a power containing [joJ • For example, we write 

g-boj-l ^ g{io}g-jo-l 

and expand separately. We then obtain the result after a long series of 
calculations and simplifications. □ 


We turn to the asymptotics of the put price. We have to replace the 
running minimum by the running maximum 


Mt := max St *, 
t*<t 


and the initial level becomes 

Jo = 


log(Mb5'i) 


(> 0 ). 


a^jrjn 
















24 


KARL GROSSE-ERDMANN AND FABIEN HEUWELYCKX 


With this re-interpretation, the numbers di, d 2 , d^, d^ and keep their 
meaning. 

For r > 0, Goldman, Sosin and Gatto m found that 

Pssit) = -Pt + StOiBi + MtB2 02)^3 

with 

Bi = ^>(di), 

B 2 = e-^^^-d 2 ), 

^3 = e— 

and, for r = 0, Babbs [1] obtained that 

PBsit) = -Pt + StBi + MtB2 + St{Bl + Bl) 


with 

B. 


3 = ( log 


Mt 


+ 




-d\/i 


Bl - (y\fT^ 


Theorem 4.3. Let 0 < t < T and n G N. The price at time t of the Euro¬ 
pean lookback put option with floating strike in the n-period CRR binomial 
model satisfies the following: 

(i) if r > 0 then 


Pl\t) = - St^ie^B^ + 02B:f)^ 




2 

St'^ ((01 + 2)5i + (02 + 2 - TOBs) + MtT2Bi 


+ 0 


where B 4 = ay/r {^) 
and T 2 = ^-\-Kn+ log 
(ii) if r = 0 then 


12 
1 

n3/2, 

^Nj(l-2r/cr2)/2 

y/ 2 'K 


1 

J n 


,T, = ^e2Kn-il + K)log 


Mt 


Mt’ 


Pl\t) = pPgit) - St^{2B, + B*, + Bl)^ 


+ 




+ 0 


a‘^T 

1 


cr^r. 


3 + - —)Bi + ^3*] + StTfBl 


1 

J n 


n 


3/2 


where Tf = + Kn + 


12 


i^log A 


dcry/r 


Mt- 


The proof is quite similar to the proof of Theorem 
omitted. 


and is therefore 


Remark 4.4. In case the option is valued at emission we have that St = Mt 
and thus Kn = 0. The formulas obtained here then reduce to the formulas 
previously found by the second author m- 
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Remark 4.5. In Theorems 14.21 and 14.31 we have that the coefficients of - 
are bounded functions of n. In fact, these coefficients are affine functions 
of Kn = {jo}(l ~ {jo})) which is bounded in n. The function x i-A x(l — x) 
is therefore responsible for the parabola-like oscillations in Figure STJ 

Remark 4.6. The coefficients for r = 0 are the limits of those for r > 0. 


5. Numerical examples 


In this part we give a numerical illustration of Theorems 14.21 and Theo¬ 
rem 1131 These results tell us that 


cl' = — + o(^) 

Jn n Xn^/^j 


1 


and 


^^(713/2) 


Pl> = pg + ^ + ^ 

n 


with certain constants Ci, Pi and functions C 2 , P 2 that are bounded in n. 
For example, for the call we should find that {Cn — and {Cn — 

Cbs ~ C\l^/n)n almost coincide respectively with Ci and C 2 for large n. 

This will be considered in the four tables below. We choose as values 
Sq = 80, fj = 0.2 and r = 1.27. For each type of option we produce an 
example with a positive value for the spot rate (r = 0.08) and an example 
with r = 0. 

For the call, we take Mt = 60 as the minimal price of the underlying (see 
Tables [T] and [5|) . The results obtained are consistent with Theorem 14.21 


Number of periods n 

1,000 

5,000 

10,000 

50,000 

100,000 



26.3647 

26.3765 

26.3794 

26.3832 

26.3842 

'-'BS 


26.3864 

26.3864 

26.3864 

26.3864 

26.3864 

{cll 

Cbs)V^ 

-0.6866 

-0.6987 

-0.7004 

-0.7040 

-0.7050 

Cl 


-0.7071 

-0.7071 

-0.7071 

-0.7071 

-0.7071 

{cl! 

~ Cbs ~ C\ly/n)n 

0.6491 

0.5931 

0.6658 

0.6868 

0.6746 

C2 


0.6640 

0.5961 

0.6635 

0.6808 

0.6681 


Table 1 

. Example for the call (r > 0) 



We take the same parameters for the put but here we consider the maximal 
price of the underlying at time t as Mt = 100 (see Tables [3] and H]) . The 
results are consistent with Theorem 14.31 

6. Conclusion 

The main goal of our paper is to derive an asymptotic expansion in powers 
of for the price of European lookback options with floating strike. In 

order to achieve this we had to refine a discrete model of Cheuk and Vorst. 
Their tree only worked for options evaluated at emission. In our work we 
consider the price of the option at any time between emission and maturity. 
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Number of periods n 

1,000 

5,000 

10,000 

50,000 

100,000 

rfi 


21.3779 

21.4016 

21.4074 

21.4151 

21.4169 

/^fi 

'^BS 


21.4214 

21.4214 

21.4214 

21.4214 

21.4214 

{ci! 


- 1.3755 

- 1.3956 

- 1.3985 

- 1.4044 

- 1.4060 

Cl 


- 1.4095 

- 1.4095 

- 1.4095 

- 1.4095 

- 1.4095 

{Cl^ 

- ^BS - Ci/^yn)n 

1.0746 

0.9868 

1.1024 

1.1371 

1.1173 

C2 


1.1144 

1.0069 

1.1136 

1.1410 

1.1209 


Table 2 

. Example for the call (r = 

0 ) 


Number of periods n 

1,000 

5,000 

10,000 

50,000 

100,000 

pf^ 

n 


16.3662 

16.4536 

16.4747 

16.5031 

16.5098 

pfi 

^BS 


16.5260 

16.5260 

16.5260 

16.5260 

16.5260 

{pI^ 

Pbs)V^ 

- 5.0523 

- 5.1200 

- 5.1274 

- 5.1325 

- 5.1394 

Pi 


- 5.1466 

- 5.1466 

- 5.1466 

- 5.1466 

- 5.1466 

{pI^ 

- Pbs - Pil^M^ 

2.9814 

1.8813 

1.9153 

3.1439 

2.2781 

P2 


3.0671 

1.9652 

1.9524 

3.1866 

2.3146 


Table 3 

. Example for the put (r > 0) 


Number of periods n 

1,000 

5,000 

10,000 

50,000 

100,000 

pf^ 

n 


23.4800 

23.5410 

23.5559 

23.5759 

23.5806 

pfi 

^BS 


23.5921 

23.5921 

23.5921 

23.5921 

23.5921 

{pI^ 


- 3.5462 

- 3.6140 

- 3.6217 

- 3.6271 

- 3.6340 

Pi 


- 3.6413 

- 3.6413 

- 3.6413 

- 3.6413 

- 3.6413 

{pI^ 

- p^^s - Pi/V^h 

3.0079 

1.9330 

1.9554 

3.1721 

2.3143 

P2 


3.0623 

1.9703 

1.9577 

3.1807 

2.3166 


Table 4. Example for the put (r = 0) 


This more general situation turned ont to require a new type of tree that 
mixes a partial binomial tree with a Cheuk-Vorst tree. Counting the number 
of paths in this tree, we derive a closed formula for the option price. 

In order to describe the asymptotic behavionr of this formula we need 
an asymptotic expansion of the binomial cumulative distribution function 
with a smaller error term than known so far in the literature. Following the 
work of Chang, Lin and Palmer i, na, we base our work on an integral 
representation of the binomial cumulative distribution function that is due 
to Uspensky m- 
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We procure explicit formulas for the coefficients of and n~^ in the 

asymptotic expansion, both for the call and the put, and for any values 
of the parameters; in particular, we allow the spot rate to be zero. These 
formulas confirm the convergence to the Black Scholes prices. Our results 
are tested on random examples. 

Several issues can be proposed as a follow-up to our work. One can 
continue the study on lookback options; so far no asymptotic expansions are 
known for the fixed-strike case. Cheuk and Vorst [1] built a one-state tree 
which can be the basis for this work. However, the price deduced from their 
tree cannot, a priori, be written using a binomial cumulative distribution 
function. It will be necessary to provide an asymptotic expansion for some 
new type of functions. 

Another possibility is to look at Asian options for which the payoff is 
determined by the average value of the underlying. The average can be 
taken to be geometric or arithmetic. Again there are fixed-strike and floating 
strike options. So far, no closed form for the price is known in the case of 
arithmetic averages [5]. Thus determining the asymptotic expansion of the 
price obtained by an equivalent CRR tree could allow to find a closed-form 
solution also in the arithmetic case. 


Appendix A. Some integrals 


We evaluate here some integrals that are needed in the proof of Theorem 
Id. 11 Recall that R 2 = —\y and a = yVV. The first integral is a Fourier 
sine transform: 


1 2 2 1 

— / — siniao?) dip = — - sm{yx) dx 

271 Jo V’ -n- Jo X 

n 1 2 1 

/ e 2 * dx = $(?/) - - 

Jo 




with the standard normal cumulative distribution function; see m pp. 
128-129], [g p. 73]. 

The remaining integrals are Fourier sine or Fourier cosine transforms of 

1 2 

the functions x i—>■ x"*e“ 2 * , m > 1. It is well known that their values 
involve the Hermite polynomials Hm'- 


1 7 “ 

TT Jo 


_ _i „2 sin(ux) if m is odd 
X e ^ ' 


cos{yx) if m is even 
see [HI p. 15, p. 74]. For m = 1 we thus have 


dx = (—1) 


[ m / 2 j 


e 2 


l ..,2 


-iv 






1 r 

Jo 


e^^^^Jipdip = 


1 1 


24F vr Jq 
1 


xe 2 ^ sin(yx) dx 


24V ^ 
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For m = 2, 

^ roo 


_L / T ,.2 J,. P-Q ^ 


27r 


'0 


J2ip^ dip = - 


eW vr Jo 


x^e 2^ cos(yx)dx 


- 1 ), 


Q\/V 

For the following values of m we obtain 


1 


dip = 


JS 1 


2F2 


TT 


roo ^ ^ 

/ x^e“2^ sin(yx) dx 

Jo 


J 3 


2F2 ^ 


y(3-y^) 


^ith J 3 * = ^ + 


7 «/4*(p-9) 1 


— J^if^ dip = 

./n 


2F3/2 


1 P 4 -1x2 7 ^ , 

— X e 2 cos(yxjdx 
TT Jo 


Jj {P - q) e 2^^ 

2F3/2 ^ 


(3 - 62/^ + /) 


with + ^(1 - 12p; 


roo 

Ytt 


_L / pR2‘P^ T^.O A,„ — *^5 


dv9 = 


J? 1 


/o 


2F2 


TT 


J 5 e-""' 


x^e 2^ sin(2/x) dx 


y(15 - lOy^ +/) 


2F2 ^ 

with J5* = ^(1 -pg) - i(^ - Ipq+p^q^) - ^V{p-qf; 


1 


Joip^ dp = 


Jpip-Q) 1 r.6„-lx2 

21J3/2 TT 

Joip-Q) 


/‘OO ^ 2 

/ x®e“2® cos{yx) dx 

Jo 


2F3/2 ^ 


(15 - 45y2 + 15/ - y 


with Jg* = -12(6-1^9); 

/*oo 




J4 1 


27r 


2F2 vr 
7 


J4 


x'e 2* sin(yx)dx 


y(105 - 105/ + 21/ - / 


21J2 ^ 

with J* = ^(1 - pg)2 + _ qf{i - I2pq) _ _ g)2; 

^ f Js/ dv9 = ^ [ x®e"/" cos(yx) dx 


1296F3/2 TT 7o 
{p — q)^ e“2?< 


12961J3/2 ^ 


i^8(/ 
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with Hs{y) = 105—420y^+210i/^ —28y®+i/®; 


2i7„ 


Jgif^ d(p = x^e 2^%in(yx)dx 


2F2 

■mp-Q? 


2F2 ^ 


-f^9(y) 


with Jg = 


1 


-Tl 4 (i -w) and i? 9 (y) = y(945-1260y2+378/-36/+/); 

i” i sin(,x) d. 


(p_g)4 

31104F2 ^ 


Hii{y) 


with Hii{y) = y(-10395+17325y2-6930/+990/-55y®+y^°). 
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